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ABSTRACT 

Wc present first results from a Keck/DEIMOS spectroscopic survey of red 
giant brancli (RGB) stars in M33. The radial velocity distributions of the stars 
in our fields are well described by three Gaussian components, corresponding to a 
candidate halo component with an uncorrected radial velocity dispersion of o" ~ 
50kms^^, a candidate disk component with a dispersion a ~ 16kms~^, and a 
third component offset from the disk by ~ 50kms~^, but for which the dispersion 
is not well constrained. By comparing our data to a model of M33 based upon 
its HI rotation curve, we find that the stellar disk is offset in velocity by ~ 
25kms~^ from the HI disk, consistent with the warping which exists between 
these components. The spectroscopic metallicity of the halo component is [Fc/H] 
~ —1.5, significantly more metal-poor than the implied metallicity of the disk 
population ([Fe/H] ~ —0.9), which also has a broader colour dispersion than the 
halo population. These data represent the first detections of individual stars in 
the halo of M33 and, despite being ~ 10 times less massive than M31 or the 
Milky Way, all three of these disk galaxies have stellar halo components with a 
similar metallicity. The color distribution of the third component is different to 
the disk and the halo, but is similar to that expected for a single, coeval, stellar 
population, and could represent a stellar stream. More observations are required 
to determine the true nature of this intriguing third kinematic component in 
M33. 

Subject headings: galaxies: individual (M33/Triangulum) — Local Group — 
galaxies: stellar content — galaxies: structure 

1. Introduction 

M33, the Triangulum Galaxy (1^33™5P, 30°39'36"), is the third brightest Local Group 
galaxy, with an integrated luminosity of My = —18.9. It is a late- type spiral, Sc II-III, with 
two open spiral arms and no clear evidence of any bulge component (Bothun 1992; Minniti 
et al. 1993; McLean & Liu 1996). The overwhelming majority of hght in M33 is distributed 
in an exponential disk component (de Vaucouleurs 1959), and the optical disk is tilted at 
nearly 30° to the strongly warped HI envelope (Rogstad et al. 1976). The mass-to-light ratio 
of its nucleus is small (M/L < 0.4), ruling out the presence of a supermassive black hole 
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(Kormendy & McClure 1993; Lauer et al. 1998) like that implied to exist in the Milky Way 
(MW; Genzel et al. 1997) and M31 (Dressier & Richstone 1988; Kormendy 1988). 

Chandar et al. (2002) analyse the kinematics of stars clusters in M33, and find two dis- 
tinct kinematic populations, which they argue represent a disk and halo population, although 
individual stars belonging to a halo component in M33 had not been directly observed. In 
a companion paper by A. Ferguson et al. in preparation (see also Ferguson et al. 2006), 
we present our wide field photometry of this galaxy taken with the Isaac Newton Telescope 
Wide Field Camera (INT WFC), and show that a low level stellar component dominates 
the radial profile of M33 at large radius, implying the presence of an extended stellar halo 
in M33. Unlike its massive neighbor, M31, there is no evidence of substructure in M33 to a 
surface brightness threshold of /ly — 30magssq.arcsec. 

In this Letter, we present first results from a spectroscopic survey of red giant branch 
(RGB) stars in M33 using the DEep Imaging Multi-Object Spectrograph (DEIMOS) on 
Keck II (Faber et al. 2003). Section 2 describes our observations. Section 3 analyzes the 
radial velocity distributions in our fields. Section 4 analyzes the metallicities of our target 
stars, and Section 5 discusses our results and concludes. We assume a distance to M33 of 
809± 24kpc (McConnachie et al. 2004a, 2005). 

2. Observations 

Since 2002, our group has been conducting a spectroscopic survey of M31 using DEIMOS 
on Keck II (Ibata et al. 2004, 2005; McConnachie et al. 2004b; Chapman et al. 2005, 2006). 
In September - October 2005, we extended the scope of this survey to M33. The two fields 
analyzed here were taken at the end of the night at high airmass, during imperfect weather 
conditions. Their location relative to the disk of M33 can be seen in Figure 1, and are 
positioned ~ 38' (9kpc) along the southern major axis to study the kinematics of stars at 
the edge of the M33 disk. In the same way as our M31 study, we used the 'mini-slitlct' 
DEIMOS setup for these fields, which significantly increases the multipexing capabilites of 
this instrument (Ibata et al. 2005; Chapman et al. 2006). Targets were selected from the INT 
survey using a color - magnitude box defined by 20.5 < i < 22.0 and 1.0 < {V — i)^ < 4.0. 
Other targets brighter than i — 22 were automatically chosen to fill in the available space of 
the spectrograph detector. 

The contours in Figure 1 show the expected mean stellar heliocentric radial velocities 
for a kinematic model of M33, based upon the HI rotation curve of M33 derived by Corbelli 
(2003) and shown in her Figure 5. We assume a stellar distribution in M33 described by 
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p{R,z) oc exp{-R/Re)exp{-z/ze), with Re = 5.8' (1.4 kpc; Regan & Vogel 1994) and 
Ze — 0.2 kpc. Changing z^ within the range [0, 1] kpc makes no significant difference to any 
of our results. Stars occupy circular orbits with Vz = and defined by interpolation of 
the HI rotation curve, linearly extrapolated beyond R = 77' (18 kpc) where necessary. The 
model is then corrected for the position angle {6 ~ 23°) and inclination {i 54°) of the 
stellar disk of M33. The mean velocity along any sight-line through the disk is calculated 
as vi = J V • L p{R, z) dlj J p{R, z) dl, and corrected to the heliocentric frame. L is a unit 
vector pointing along the direction of our line of sight. 

The data were processed and reduced in the same way as the main M31 survey (Ibata 
et al. 2005; Chapman et al. 2006). Due to the generally poor observing conditions, the data 
analyzed here consist of 280 stars with a cross correlation coefficient > 0.02, a signal-to-noise 
(S/N) in the continuum > 2, a Tonry-Davis coefficient (Tonry & Davis 1979) > 2, and a 
velocity error Verr < 25kms~''^ {verr — lOkms"-^). 

3. Kinematics 

The top two panels of Figure 2 shows the hehocentric radial velocity (vh) distributions 
in our two fields. The dashed lines indicate Vh — 0, and the dot-dashed fines indicate the 
systemic velocity of M33, Vh — — 179kms~^. The dot-dashed histograms show the expected 
hehocentric velocity distributions of the MW stars that lie in each field using the MW model 
of Robin et al. (2003), and which satisfy our selection criteria for M33 RGB stars. The 
histograms have been scaled to match the area of the DEIMOS fields. They show that any 
contribution from MW foreground stars at Vh < — 80kms~^ is small for the purposes of this 
Letter, and so this velocity cut is adopted in the following kinematic analysis. Stars with 
Vh > — 80kms~^ are predominantly MW foreground, and the small numbers we observe 
further suggest that the contamination at more negative velocity will be negligible. 

The third panel of Figure 2 shows the velocity distributions of Fields 157 and 158 
added together to boost statistics. We expect the main stellar components in M33 will 
be well fit by individual Gaussians, and any stellar halo component in M33 will have a 
velocity distribution centered on Vh — — 179kms~^. Since M33 is ~ 10 times less massive 
than M31, we expect that the M33 halo velocity dispersion will be ~ VTO that of M31 
(cm3i ~ 145kms~^; Chapman et al. 2006). The broad Gaussian overlaid in the third panel 
of Figure 2 represents a M33 halo component with a dispersion of cr = SOkms"^, scaled to 
match the approximate number of stars observed. We also expect that the prominent peak 
of stars in these fields is due to the dominant disk of M33. Therefore, we have subtracted the 
halo Gaussian from the Vh distribution in the third panel and fitted a Gaussian to this peak. 
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centered at vi — — 102±2 km with an uncorrected dispersion ai — 16±3 km s~^ (corrected 
dispersion of ~ 12.5 kms~^). The residuals of this fit show a prominent grouping of stars at 
V ~ — 150kms~^. We fit a third Gaussian to this feature, with a peak V2 — — 150±4kms~^ 
and a dispersion cr2 = 10 ± 4kms~^. The latter quantity is more poorly constrained than 
given by the formal error, however, as its dispersion is covariant with the dispersion of the 
disk and halo, which have been fitted independently. All 3 Gaussians and their summed 
total are overlaid in the third panel of Figure 2. 

The fourth panel of Figure 2 shows the expected distribution of stars at the positions 
of Fields 157 and 158, for the model of M33 shown in Figure 1, combined with a halo model 
with a radial velocity dispersion of 50kms^^. The radial velocity dispersion of the disk is set 
to 16kms~^, to mimic our observations. The normalization of both components is chosen to 
approximately match the corresponding peaks in the data . A comparison between model 
and data shows (1) the peak at V2 — — 150kms~^ does not naturally arise as a result of our 
line-of-sight through the main M33 disk sampling stars at different circular velocities, (2) 
the stellar disk is offset in Vh from the HI disk by ~ 20 — 25kms~^, implying a dcprojected 
offset of ~ 25 - 30kms-^ For comparison, the asymmetric drift between the stars and HI 
in the MW and M31 is of order 15 — 20kms~^. Scaling for mass, the expected asymmetric 
drift in M33 is not likely to be more than ~ lOkms"^. We believe that the velocity offset 
between the two disks is due to differential warping, specifically that the HI disk is warping 
into our hne-of-sight at large radius from M33. This has been implied, quite spectacularly, 
by Rogstad et al. (1976) (see their Figure 14), and is consistent with the magnitude of the 
velocity offset. 



4. Metallicities 

Figure 3 shows the color - magnitude diagrams (CMDs) for the spectroscopically ob- 
served M33 RGB stars in Fields 157 and 158, divided by the radial velocity cuts indicated 
so as to kinematically separate the three candidate components as best as possible. Victoria 
- Regina isochrones from VandenBerg et al. (2005), with BVRI color - Tg// relations as 
described by VandenBerg & Clem (2003), corresponding to a 13 Gyr stellar population with 
metaUicities of [Fe/H] = —2.0, —1.3, —0.8 and —0.5, are also overlaid. 

The differences between the CMDs in Figure 3 are striking, and appear to validate our 
kinematic deconvolution. The candidate M33 halo stars are predominantly bluer and have 
a smaller color dispersion than the candidate disk stars, indicating a difference in the age 
and/or metallicity properties of the populations. Assuming both populations arc 13 Gyr s 
old implies a median photometric metallicity for the halo of [Fe/H] = — 1.3 ± 0.1, with 
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a dispersion of (7[Fe/H] = 0.21. Given the small number of stars which contribute to this 
measurement, this is consistent with the Ferguson et al. (2006) metallicity measurement for 
the halo, of [Fe/H] = —1.5 ± 0.1. The disk has an implied photometric metallicity of [Fe/H] 
= —0.9 ±0.1 and a broader dispersion of o"[Fe/H] = 0.35. The median color (metallicity) 
of the 'third component' is similar to the disk, although its color dispersion is remarkably 
narrow (crpe/H] = 0.15). This population is difficult to isolate kinematically, and our sample 
could be contaminated by both of the other two components. However, the locus of its stars 
on the CMD does not clearly associate it with either the candidate disk or halo population, 
and is reminiscent of a CMD for a single stellar population. If this population was due to 
MW foreground contamination, we would not expect its stars to occupy so well defined a 
locus in this diagram. The most probable explanation is that these stars belong to a stellar 
stream. These stars do not, however, define a spatially obvious structure, but given our 
limited field size and that our choice of slitlet position will affect the spatial distribution, 
this test is inconclusive at best. 

We calculate the spectroscopic metallicity of the halo stars by measuring the equivalent 
width of the CaT from stacked spectra, using an identical method to Ibata et al. (2005) 
and Chapman et al. (2006). We find that [Fe/H] = —1.5 ± 0.3, in good agreement with 
the photometric measurement. The large uncertainty refiects the relatively low signal-to- 
noise of the stacked spectrum. We are unable to measure the metallicities of the other two 
components in this way, since their velocities cause the CaT to lie on top of strong sky lines 
which prevents accurate equivalent width measurements. 

5. Discussion 

M33 is considered an example of the quintessential disk galaxy, with no evidence for 
any bulge or spheroid component. However, the radial profile for M33 presented in Ferguson 
et al. (2006), based upon our INT-WFC observations, reveals for the first time the presence 
of a spatially extended, low density component. Using Keck/DEIMOS, we have identified 
a population of RGB stars which have velocities centered on the systemic velocity of M33 
and which appear to have a large dispersion {a ~ 50kms~^). These stars are bluer than 
other RGB stars in M33, with a small color dispersion. Their average metallicity is [Fe/H] 
~ —1.5, and we attribute them to the newly discovered halo component of M33. The haloes 
of the MW and M31 have similar metallicities to this (Eggen et al. 1962; Chapman et al. 
2006; Kalirai et al. 2006), even although both these galaxies are ~ 10 times more massive 
than M33. It is important to emphasise that there is no evidence of a trend linking halo 
metallicity to galaxy mass or luminosity for the Local group galaxies, as has been suggested 
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from recent Hubble Space Telescope studies of galaxies within ~ 10 Mpc (Mouhcine et al. 
2005a,b). Instead, all the large Local Group spiral galaxies appear to have very similar halo 
metallicities. 

We also observe stars in the disk of M33. Their Vh distribution is well described by 
a Gaussian with an uncorrected dispersion of o" ~ 16kms~^, corresponding to a corrected 
dispersion of ~ 12.5 km s~^. These stars are redder and have a broader color dispersion 
than the halo stars. The stellar disk is offset in Vh from the HI disk in our simple model by 
~ 20 — 25kms^^, which is fully consistent with the extreme warping of the HI disk implied 
by Rogstad et al. (1976). 

Finally, we present evidence for a third kinematic feature in M33, revealed by the Vh 
distribution of the southern fields. The velocity dispersion of this component is not well 
constrained by our data. A comparison to our simple model of the disk of M33 shows that 
this feature does not arise as a result of our line-of-sight through M33 sampling stars at 
different circular velocities in the main disk. The locus of these stars in the CMD has a 
remarkably small color dispersion, reminiscent of a single stellar population and inconsistent 
with belonging to a foreground population. We believe the most likely explanation for this 
feature is that it is a stellar stream. However, more spectroscopic observations of fields 
at different locations around the disk of M33 are required before the true nature of this 
intriguing feature is revealed. 

The data presented herein were obtained at the W.M. Keck Observatory, which is oper- 
ated as a scientific partnership among the California Institute of Technology, the University 
of California and the National Aeronautics and Space Administration. The Observatory 
was made possible by the generous financial support of the W.M. Keck Foundation. We 
thank E. Corbelli for supplying us with her data for the HI rotation curve of M33, and the 
anonymous referee for constructuve comments . AWM would like to thank J. Navarro and 
S. EUison for financial support. AMNF is supported by a Marie Curie Excellence Grant 
from the European Commission under contract MCEXT-CT-2005-025869. GFL acknowl- 
edges support through ARC DP0343508. 
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Fig. 1. — A tangent plane projection of RGB stars in M33 with I < 22, from our INT WFC 
survey (Ferguson et al. 2006). The locations of the Keck/DEIMOS fields are indicated. 
Contours show the mean hehocentric radial velocity of a model M33 disk, based upon the 
CorbeUi (2003) HI rotation curve. 
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Fig. 2. — First and second panels: heliocentric radial velocity distributions for our M33 fields. 
The dashed lines represent zero velocity, and the dot-dashed line represents the systemic 
velocity of M33, Vh = — 179kms~^. The dot-dashed histograms show the heliocentric velocity 
distribution of the expected MW foreground stars in each field, from the MW model by Robin 
et al. 2003, which satisfy our color selection criteria. These histograms have been scaled to 
match the area of the DEIMOS fields, and the small number of stars at Vh > — 80kms^^, 
which are most likely foreground, suggest that contamination at more negative velocity 
will be negligible. Third panel: Same as first two panels, except now the two fields have 
been added together to boost statistics. Three Gaussian velocity distributions and their 
summed total are overlaid. See text for details. Fourth panel: the expected heliocentric 
radial velocity distribution of stars in our fields, using our M33 disk model and an additional 
Gaussian component, representing a smooth stellar halo. 
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Fig. 3. — CMDs of spectroscopically observed RGB stars from Fields 157 and 158, selected to 
satisfy the velocity cuts indicated (ai = 16kms~^, 02 = lOkms"^). These velocity cuts pre- 
dominantly sample stars which belong to the candidate halo (left panel), disk (middle panel) 
and the 'third component' (right panel). Victoria - Regina isochrones (VandenBerg et al. 
2005), corresponding to 13 Gyr old stars and with metallicities of [Fe/H] = —2.0, —1.3, —0.8 
and —0.5, are overlaid. 



